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We report the first reconstruction in hadron collisions of the suppressed decays B~ — > D{— > 
K + n~)K~ and B~ — > D(—¥ K + n~)n~ , sensitive to the CKM phase 7, using data from 7 fb _1 of 
integrated luminosity collected by the CDF II detector at the Tevatron collider. We reconstruct a 
signal for the B~ — > D(— > K + tv~)K~ suppressed mode with a significance of 3.2 standard deviations, 
and measure the ratios of the suppressed to favored branching fractions R(K) — [22.0 ± 8.6(stat) ± 
2.6(syst)] x 10~ 3 , R + {K) = [42.6 ± 13.7(stat) ± 2.8(syst)] x 10~ 3 , R~ (K) = [3.8 ± 10.3(stat) ± 
2.7(syst)] x 10~ 3 as well as the direct CP-violating asymmetry A(K) = -0.82±0.44(stat)±0.09(syst) 
of this mode. Corresponding quantities for B~ — > D{—> K + ir~)ir~ decay are also reported. 

PACS numbers: 13.25.Hw 11.30.Er 14.40.Nd 



'Deceased 

t\Vith visitors from "Istituto Nazionale di Fisica Nucleare, Sezione 
di Cagliari, 09042 Monserrato (Cagliari), Italy, ^University of CA 
Irvine, Irvine, CA 92697, USA, c University of CA Santa Barbara, 
Santa Barbara, CA 93106, USA, ^University of CA Santa Cruz, 
Santa Cruz, CA 95064, USA, e CERN,CH-1211 Geneva, Switzer- 
land, -f Cornell University, Ithaca, NY 14853, USA, ^University of 
Cyprus, Nicosia CY-1678, Cyprus, ^Office of Science, U.S. Depart- 
ment of Energy, Washington, DC 20585, USA, "University College 
Dublin, Dublin 4, Ireland, -'University of Fukui, Fukui City, Fukui 
Prefecture, Japan 910-0017, k Universidad Iberoamericana, Mex- 
ico D.F., Mexico, ; Iowa State University, Ames, IA 50011, USA, 
""University of Iowa, Iowa City, IA 52242, USA, "Kinki University, 
Higashi-Osaka City, Japan 577-8502, "Kansas State University, 
Manhattan, KS 66506, USA, ^University of Manchester, Manch- 
ester M13 9PL, United Kingdom, 'Queen Mary, University of Lon- 
don, London, El 4NS, United Kingdom, ""University of Melbourne, 
Victoria 3010, Australia, s Muons, Inc., Batavia, IL 60510, USA, 
'Nagasaki Institute of Applied Science, Nagasaki, Japan, "National 



The measurement of CP-violating asymmetries and 
branching ratios of B~ —> DK~ [l[ decay modes 
allows a theoretically clean extraction of the phase 
7 = &rg(~V ud V* b /V cd V* b ) of the Cabibbo-Kobayashi- 
Maskawa quark- mixing matrix Vqkm , & fundamental pa- 
rameter of the standard model [4]. In these decays the 
interference between the first order tree amplitudes of 
the b — » cus and b — > ucs processes leads to observables 
that depend on their relative weak phase 7, their relative 
strong phase 8b, and the magnitude of the amplitude ra- 
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tio tb 9 ■ These quantities can all be extracted from data 
by combining several experimental observables. This can 
be achieved in several ways, using a variety of D decay 
channels [H-IE]- An accurate knowledge of the value of 
7 is instrumental in establishing the possible presence 
of additional non-standard model CP-violating phases in 
processes where higher-order diagrams are involved 0, Hi- 
lts current determination has a relative uncertainty, dom- 
inated by statistical uncertainties, between 15 and 20%, 
depending on the method [||. 

A promising class of processes consists of B meson de- 
cays that are a coherent superposition of the color fa- 
vored B~ — > D°K~ followed by the doubly Cabibbo sup- 
pressed d_ecay D° — » K + n~ , and of the color suppressed 
B~ — > D°K~ followed by the Cabibbo favored decay 
D° — > K + tt~. The magnitude of the two amplitudes is 
comparable, allowing for large CP-violating asymmetries 
sensitive to_the phase 7. The following observables can 
be defined [f|: 



R(K) 



R±(K) 



A(K) 



B{B~ 



[K+n-] D K-)+B{B+ -> [K-tt+] d K^ 



B(B~ 
B(B ± 



[K-w+] D K-) + B{B+ 
[K^Tr±] D K±) 



[K+ir-] D K+)' 



B(B ± 
B(B~ 



[K+tt-] d K-) - 



B(B+ 



\K- 



~}dK+) 



where B 
mode and 



B(B~ -> [K+ir-] D K-) + B(B+ -> [K~it+] d K+) ' 

— > [K + ir~]]jK~ is the suppressed (sup) 
B~ — s- [K~tt + ]£)K~ is the favored (fav) 
mode. In the approximation of negligible CP-violation 
in D decays and negligible D° — D° mixing, whose ef- 
fects were shown to be small in Ref. [13], these quan- 
tities are related to the CKM phase 7 by the equa- 
tions 

R ± : 



1 R = 



1 D ' 1 B ~ 

^■tbtd sin 7 sin (8b 



Sd)/R, where 



d + r B + 2r D r B cos7Cos(<5 B + S D ), 
- 2tdtb cos (5b + So ± 7), and A = 

A(D"^K + 7r- 
A(D u ^K-ir+ 

and 80 is the corresponding relative strong phase. The 
smallness of the product of branching fractions for these 
suppressed final states (0(1O -7 )) has been a strong limi- 
tation to their use in 7 determinations. Evidence for the 
suppressed B~ — > DK~ channel has only recently been 
obtained by the Belle collaboration [ll[. The large pro- 
duction rate of B mesons available at hadron colliders of- 
fers a unique opportunity for improving the experimental 
determination of the angle 7. Measurements of branching 
fractions and CP-violating asymmetries of B~ — > DK~ 
modes in less suppressed final states of the D meson ( CP- 
even modes K~K + and tt~tt + ) have already been per- 
formed in hadron collisions [l2[. However, the small de- 
cay rates along with large potential backgrounds from 
misidentified favored decays, which only differ for the 
identity of the final particles, make the reconstruction of 
suppressed modes in hadron collisions significantly more 
challenging. 

In this Letter, we describe the first reconstruction of 
B~ — > D sup K~ modes performed in hadron collisions, 
based on data from a total integrated luminosity of 7 fb _1 
of pp collisions at = 1.96 TeV, collected by the up- 
graded Collider Detector (CDF II) at the Fcrmilab Teva- 
tron. We report measurements of R(K), R ± (K), and 



A(K) for those modes. We also report measurements re- 
lated to the corresponding Dn~ modes, since measurable, 
albeit smaller, 7-dependent asymmetries may also be 
found in these modes Q . The maximum possible value of 
the asymmetry is A max = 2rBrD/(r%+rf :) ), where tb can 
be rs(K) or rs(7r). Taking into account the CKM struc- 
ture of the contributing processes, we expect that rsi^) 
is suppressed by a factor |V r c( jV r qts / 1 K l( iV r cs | ~ tan 2 0q with 
respect to tb(K), where 8c is the Cabibbo angle, and 
we assume the same color suppression factor for both 
DK and Dtt modes. Using r B (K) = 0.103±^ 0], 
r B (7r) - 0.005 @, and r 2 D = (3.80±0.18) x 10" 3 0], we 
expect A max (K) rj 0.90 and A max (n) w 0.16. 

CDF II is a multipurpose magnetic spectrometer sur- 
rounded by calorimeters and muon detectors, and is de- 
scribed in detail elsewhere [l4| 03 • The resolution on 
transverse momentum of charged particles is u Pt /pt — 
0.07%pT/(GeV/c), corresponding to a typical mass res- 
olution of 18 MeV/c 2 for our signals. The specific ion- 
ization energy loss dE / dx of charged particles can be 
measured from the charge collected by a gaseous drift 
chamber (COT), and provides 1.5er separation between 
pion and kaon particles for p > 2 GeV/c. Candidate 
events for this analysis are selected by a three-level on- 
line event-selection system (trigger). At level 1, charged 
particles are reconstructed in the COT b y a hardware 
processor, the extremely fast tracker XFT Two op- 
positely charged particles are required, with transverse 
momenta pr > 2 GeV/c and scalar sum pxi +PT2 > 5.5 
GeV/c. At level 2, another processor, a silicon vertex 
trigger (SVT) fl9| . associates r — <j> position measure- 
ments from an inner silicon detector with XFT tracks. 
This provides a precise measurement of the track impact 
parameter do, the transverse distance of closest approach 
to the beam line. The resolution of the impact parame- 
ter measurement is 50 ^m for particles with px of about 
2 GeV/c, including a sa 30 (im contribution due to the 
transverse beam size, and improves for higher transverse 
momenta. 

We select B hadron candidates by requiring two SVT 
tracks with 120 < do < 1000 /urn. To reduce background 
from light-quark jet pairs, the two trigger tracks are re- 
quired to have an opening angle in the transverse plane 
2° < Atp < 90°, and to satisfy the requirement L xy > 200 
/_tm, where L xy is defined as the distance in the transverse 
plane from the beam line to the reconstructed two-track 
vertex. The level 1 and 2 trigger requirements are then 
confirmed at trigger level 3, where the event is fully re- 
constructed in software. 

The events collected by the trigger are further selected 
by searching for a pair of oppositely charged particles 
compatible with a two-body D decay. The invariant mass 
AId of the pair is reconstructed for both pion and kaon 
assignments of particle identities. Events are accepted for 
the analysis only when one of the possible masses is com- 
patible with the nominal D mass 1.8495 < M D < 1.8815 
GeV/c 2 , and the alternative combination, Msw(D), is 
outside a veto region of 1.8245 < M S w(D) < 1-9045 
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GeV/c 2 around the nominal D mass. The D candidate 
is then combined with a negatively charged particle in the 
event with pt > 0.4 GeV/c to form a B~ candidate. A 
three-dimensional kinematic fit of each decay candidate 
trajectory is performed by constraining the two tracks 
forming the D candidate to a common vertex and to the 
nominal D mass; the D candidate and the remaining 
track to a separate vertex; and the reconstructed momen- 
tum of the B~ candidate to point back to the primary 
pp interaction vertex determined from other tracks in the 
event. 

The events are then divided into two non-overlapping 
samples, nominally classified as favored or suppressed, 
according to the relative charge of the B candidate with 
the decay product of the D that has been classified as the 
kaon. The veto requirements applied to the D mass re- 
constructed with the alternative particle assignment re- 
move a large fraction of the background of favored de- 
cays from the sample classified as suppressed, and vice 
versa, ensuring no overlap between the samples and a 
complete symmetry of the selection, which is a crucial 
aspect of the analysis. The small residual contamination 
of each sample from events with an incorrect identifi- 
cation of D decay products is accounted for as part of 
the inclusive background B~ — > !?(—>• X)ir~ , where X 
are modes other than Kir (see below). A further veto is 
applied to the invariant mass formed by the track from 
the B candidate and the oppositely charged track from 
the D candidate, again requiring it to be incompatible 
with the D meson mass, using the same range as the 
first veto. This requirement suppresses the contamina- 
tion from tracks from real B decays that have been in- 
correctly labeled as D decay products, and is applied 
symmetrically to both samples. A further suppression of 
this background is achieved by requiring that the trans- 
verse distance between B and D decay vertex is greater 
than 100 fim. This has the additional effect of reducing 
contamination from non-resonant three-body decays of 
the type B + — > h + h~h + , in which all tracks come from 
a common decay vertex, and where h indicates either K 

Or 7T. 

Additional requirements are applied to the following 
observablcs: the impact parameter ds of the recon- 
structed B candidate relative to the beamline; the iso- 
lation of the B candidate I B [20[ ; the goodness of fit of 
the decay vertex Xb', the significance of the B hadron 
decay length L xy (B) /ctl ^b)', the angle a between the 
three-dimensional momentum of the B candidate and the 
three-dimensional decay length; AR = yf A<p 2 + Arj 2 be- 
tween the track from the B hadron and the D meson; 
the cosine of the angle between the D and the flight di- 
rection of the B, in the B meson rest-frame, cos (91k; the 



difference of the kaon probability [21| values of the tracks 
forming the D to discriminate kaon-pion pairs from pion- 
pion and kaon-kaon pairs, An. The threshold values 
for all these requirements, and for the allowed D mass 
window mentioned above, were determined by an unbi- 
ased optimization procedure, maximizing the quantity 



Ns/(1.5 + tJNb) [22|], with no use of simulated signal. 
The signal Ns is defined as the expected rate of sup- 
pressed B~ — > D sup TT~ events. We take advantage of our 
large sample of favored B~ — > Df av ir~ decays, using it 
as a model for the kinematical and particle identification 
properties of the suppressed decay by simply considering 
the swap in sign. The resulting requirements are the fol- 
lowing: L X y(B)/a Lxy{ B) > 12, d B < 50 /im, x% < 13 > 
J B (cone = 1) > 0.4, 7 B (conc = 0.4) > 0.7, a < 0.15, 
AR < 1.5, | cos 0£| < 0.6, Are > -1. After applying all 
the above selection criteria, the invariant mass of each 
B~ — > Dh~ candidate is evaluated using a nominal pion 
mass assignment to the particle h~ coming from the B 
decay. Figure Q] shows the distributions for B ± candi- 
dates. 

With the help of large simulated samples of B mesons, 
we determine that the only modes contributing non- 
negligible backgrounds are B~ — > D{— > X)h~, B~ — > 
D*°ir~, with D*° -> D ^/^ , non-resonant B~ -> 
K~ir + ir~, and B Q — > D*^~l + vi. The large contribution 
of B~ -> D{-^ K + I\-)h- reported in Ref. [U HI is 
strongly suppressed by our selection, since we reconstruct 
the D mass in the Kn mass hypothesis. 

We use an extended unbinned maximum likelihood fit, 
exploiting mass and particle identification (PID) infor- 
mation to statistically separate the B~ — > DK~ and 
B~ — > Dtt~ signals, the combinatorial background, and 
the physics backgrounds. PID information on the track 
from the B decay is incorporated in the kaon probabil- 
ity observable [2l|. The extended likelihood function is 
defined as C = YliViCi, where i runs over the favored 
and suppressed modes, positive and negative charges. 

N tot 

The Poisson distribution V% is equal to ^ N t t x e~ M , where 

Nf ot is the number of events of each sub-samples and 
fi is the expected mean value. The individual likeli- 
hood components have the following structure: £j = 

Fir * fjPj(M r , K r \9 r ), where / and P(M r , n r \9 r ) are 
the fractions and the probability density functions of the 
signal and background modes, and 6 r are other free pa- 
rameters of the fit, a mass scale parameter with respect 
to the nominal B-mass and a scale factor multiplying the 
width of the shapes of the B~ — > D(—> K + ir~)h~ signals. 
The fit is simultaneously performed on the favored and 
suppressed samples. Common parameters are the expo- 
nential function for the combinatorial background, whose 
normalization and slope are determined by the fit; the 
functional expression for signal and background modes; 
and the ratio between B~ — > D*°ir~ and B~ — > Dir~ 
fractions. The numbers of events and the fractions of 
signal and background arc determined by the fit and the 
observables are extracted from them. We tested on sim- 
ulation that our fit does not exhibit any significant bias. 

The shape of the mass distribution assigned to each sig- 
nal and physics background has been modeled using sim- 
ulated events including the effect of final state QED ra- 
diation and parameterized with different functions. Sys- 
tematic uncertainties are assessed by varying the values 
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FIG. 1: Invariant mass distributions of B — > Dh ± for the suppressed mode (bottom meson on the left and antibottom on the 
right). The pion mass is assigned to the charged track from the B candidate decay vertex. The projections of the common 
likelihood fit (see text) are overlaid. 



of those function parameters within their errors. 

A large sample of D* + — > D°(— > K~tt + )tt + decays is 
used to calibrate the average dE/dx response of the de- 
tector to kaons and pions, using the charge of the pion in 
the D* + decay to determine the identity of the D decay 
products. The shape of the k distribution is calibrated 
within our own sample, by using kaons and pions from 
the decay of the D meson in the favored sample. Un- 
certainties on the calibration parameters are included in 
the final systematic uncertainty of A, R and R^, taking 
into account the full correlation matrix of the par ameters 
characterizing the shape of the k distribution[23(. 

The B~ — >• DK~ and B~ -> Dit~ event yields ob- 
tained from the fit to the data are reported in Table Q] 
Fit projections on the invariant mass distributions are 
given in Fig. [T] They provide a consistent description 
of the observed distributions in the data. We find evi- 
dence for a signal in the B~ — > DK~ suppressed mode 
with a significance of 3.2a. The significance is evalu- 
ated by comparing the likelihood-ratio observed in data 
with the distribution expected in statistical trials. Sev- 
eral distributions are generated corresponding to differ- 
ent choices of systematic parameters. The quoted signif- 
icance corresponds to the distribution yielding the most 
conservative p-value. The raw fit results are then cor- 
rected for the reconstruction efficiency e, due to differ- 
ent probabilities of K + , K~, tt + and ir~ to interact 
with the tracker material. We use previous measure- 
ments of #S = 1.0178 ± 0.0023(stat) ± 0.0045(syst) 

and !0l = 0.997 ± 0.003(stat) ± 0.006(syst) 0. We 

extract = 0.998 ± 0.015(stat) ± 0.016(syst) from 

our own sample of favored B~ — > Dn~ decays. 



Systematic uncertainties are determined by repeating 
the fit changing the mass and the dE/dx model (Ta- 
ble [nj. The dominant contribution is the uncertainty on 
the B~ — > D(—> X)tt~ shape. This is the largest physics 
background, and it lies under the signal peak. 

In summary, we find evidence for the B~ — > D(—> 
K + tt~)K~ suppressed mode with a significance of 3.2 
Gaussian standard deviations. We measure the ra- 
tios of the suppressed ([K + it~]e i K~ to favored 
([K-tt+] d K- /n-) branching fractions R{K) = [22.0 ± 
8.6(stat)±2.6(syst)]xl0" 3 , R + (K) = [42.6±13.7(stat)± 
2.8(syst)]xl0- 3 , R-(K) = [3.8±10.3(stat)±2.7(syst)j X 
lO- 3 and R(tt) = [2.8 ± 0.7(stat) ± 0.4(syst)] x 10" 3 , 
R+{tt) = [2.4 ± l.O(stat) ± 0.4(syst)] x 10" 3 , R-{n) = 
[3.1 ± l.l(stat) ± 0.4(syst)] x 10" 3 as well as the direct 
CP-violating asymmetries 

A{K) = -0.82±0.44(stat)±0.09(syst), 
A(tt) = 0.13±0.25(stat)±0.02(syst). 

The observed asymmetry A(K) deviates from zero by 2.2 
standard deviations. 

These measurements, performed here for the first time 
in hadron collisions, are in agreement with previous mea- 
surements from BaBar [2(| and Belle [ll[ with compa- 
rable uncertainties. These results can be combined with 
other B~ — > DK~ measurements to improve the deter- 
mination of the CKM angle 7. 

Acknowledgments 

We thank the Fermilab staff and the technical staffs 
of the participating institutions for their vital contribu- 



7 



TABLE I: B — > DK and B — > Dn event yields obtained from the fit to the data. Only statistical uncertainties are 
quoted. 



D mode 


B+ - 


> Dtt+ 


B~ -> Dtt- 


B+ -> DK+ 


B~ -> Dif ~ 


K~ty + (favored) 


9882 


± 103 


9892 ± 103 


694 ± 39 


767 ± 41 


K + -k~ (suppressed) 


24 


±9 


31 ±10 


29 ±9 


3±8 



TABLE II: Summary of systematic uncertainties. 



Source 


R(ir) 






R(K) 


R+(K) 


R-(K) 


A{k) 


A(K) 


dE 1 dx model 


<0.0001 


<0.0001 


<0.0001 


0.0001 


0.0003 


0.0001 


<0.01 


<0.01 


B~ — > D[— > X)tt~ shape 


0.0004 


0.0004 


0.0004 


0.0025 


0.0026 


0.0026 


0.01 


0.09 


Other backgrounds 


<0.0001 


<0.0001 


<0.0001 


0.0006 


0.0006 


0.0005 


<0.01 


0.02 


Efficiency 


<0.0001 


<0.0001 


<0.0001 


0.0003 


0.0009 


0.0001 


0.01 


<0.01 


Total 


0.0004 


0.0004 


0.0004 


0.0026 


0.0028 


0.0027 


0.02 


0.09 
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